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Abstract 


Optimum design of modem ground -satellite communication 
systems requires the knowledge of rain-induced differential attenuation, 
differential phase shift, and cross-polarization factors. After a 
rather comprehensive assessment: of different available analytical 
techniques, an efficient technique is chosen to yield the 
desired scattering results. Tabulation of the scattered fields 
for Pruppacher-Pitter raindrops with sizes ranging from 0.25 mm 
to 3.5 mm at 20° C and at 30 GHz for several incidence angles is 
given. 
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I. Introduction 


The optimum design of modem ground -satellite conammlcatlon systems 

1-5 

requires detailed knowledge of the rain-induced radio propagation affects. 
Investigation by Pruppacher and Fitter^ on the shapes of realistic rain- 
drops shows that raindrops take on shapes which deviate significantly from 
spheres when their sizes exceed about 1 mm. In fact their shapes even depart 
noticeably from oblate spheroidal shapes. It Is well known according to the 
electromagnetic scattering theory that the shape of a scatterer Is a crucial 
factor In determining the resonant scattering characteristics of a particle. 
Hence, accurate treatment of the scattering of realistic non-spherlcal shape 
raindrops is essential In the determination of rain-induced radio propagation 
effects which may include the generation of cross polarized fields^ and the 
introduction of Incoherent fields due to the presence of multiple scattering 
phenomenon.^ 

Various aspects of the rain scattering problems have been treated by 

a number of investigators using various Justifiable or sometimes marginal 
2-5 

techniques. The purpose of this paper is twofold. We shall first pro- 
vide a rather comprehensive assessment of previously used techniques. Based 
on our assessment we shall then select the most attractive method for raln- 
scatterlng calculation. Then, complete tabulation of scattered fields for 
Pruppacher -Pit ter varying-shape raindrops with sizes ranging from 0.25 mm 
to 3.5 mm at 2(fi C and at 30 Olz for several incidence angles will be given. 
These extensive single particle scattering results also represent the 
necessary input for the determination of the distribution of incoherent and 

g 

coherent intensities due to multiple scattering effects. 
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II. An Atsessment cf Available Techniques fot the Raindrop Scattering 
Problems . 

The classical solutions of electromagnetic wave scattering are 
limited to simple objects of separable boundaries such as spheres, cylinders, 
etc. Recent developments In scattering analysis for scatterers of arbi- 
trary shapes Include the geometric theory of diffraction for high frequency 
scattering ^ and the method of moment, the perturbation method, the 

point matching method,^ the extended boundary condl*’.ions method,^ and the 

12 

finite element technique for resonant region scattering. 

To achieve a proper perspective of the state-of-the-art of various 
analytical /numerical techniques, an assessment Is made in the following: 

(a) Geometrical Theory of Diffraction. 

This mettiod Is mainly a high frequency technique which treats the 
smooth part of the scatterer by geometric optics and obtains the scattered 
fields from edges and glazing surfaces via solutions of canonical problems 
such as scattering by wedges, tips, etc. The geometric theory of diffraction 
may yield reasonable results even at resonant frequencies. However, it is 
still a technique most conveniently applied to convex and perfectly con- 
ducting targets at high frequencies. When the targets are concave or 
transparent, tracing the multiple scattering rays makes the method overly 
cumbersome to use. 

Advantages: 

Easy to learn 

Clear physical concept 

Disadvantages: 

DlffiC'ilt to apply to dielectric bodies. 

Questionable validity at lower frequencies. 
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(b) Ite'chod of Moment 

This method is based on the solution of an Integral equation. The 
equation Is a surface Integral equation If the scatterer Is a perfectly con- 
ducting body, a pair of couvled surface integral equation If the scatterer 
is a homogeneous dielectric body, and a volume Integral equation If the 
scatt^er Is an Inhomogeneous dielectric body. If the dielectric scatterer 
has a volume greater than (O.IX)^, the deuand of the method of moment for 
computer memory or time becomes excessive. 

Advantages: 

Efficient for thin conducting scatterer. 

Disadvantages: 

Limited to dielectric scatterer of small volume. 

(c) Perturbation Method. 

This method is an extension of the Mie scatte* Ing solution to scatterers of 
non-spherical geometry by the boundary perturbation technique. This technique 
Is inher«itly limited to particle shapes which are small perturbations of a 
sphere. 

Advantages: 

Analytic expressions are available. 

Disadvantages: 

Uncertain convergence 

Limited to bodies that are small perturbations of a sphere. 

Analytic expressions for higher-order terms are excessively 
complicated. 

(d) Point Matching Me«.ho<’. 

Field expressions for the interior region of the scatterer and for the 
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exterior region are given In terms of complete sets of spherical harmonics. 
Boundary conditions are satisfied by point -matching on the bou.nlary surface. 
Advantages: 

Conceptually simple 
Disadvantages^ 

Uncertain convergence 
Demands excessive computer time 

Limited to bodies that are close to spherical shapes 

(e) Extended Boundary Condition Method. 

Integral representations for the fields can be derived which satisfy 
the wave equation and all necessary boundary conditions. By expanding the 
fields in terms of a complete set of spherical harmonics and by making use 
of analytic continuation techniques, one may reduce the Integral representations 
to a set of linear algebraic equations. 

Advantages: 

May be used to treat scatterers of large volume and arbitrary 
shapes 

Numerical results can be generated very efficiently 

Disadvantages: 

Usually not applicable to Inhomogeneous dielectric scatterers 

(f) Global-Local Finite Element Method. 

The term Global-Local Finite Element Method refers to a numerical analysis 
technique in which both the contemporary finite element and classical Saylelgh- 
Bltz approximations are employed. The scattering problem Is divided Into two 
regions: the exterior region and the Interior region. Without any sacrifice 
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of rigor. It Is assumed that the boundary surface between the exterior and 
Interior regions Is a sphere for three •dimensional problems acd a circular 
cylinder fcr two-dimensional problems. It is further assumed that the 
boundary surface Is so chosen that the medium In the exterior region Is 
nomogeneous and that the Inhomogoneous Irregularly-shaped absorbing object 
Is contained within the Interior region. Solutions of wave equations L. the 
exterior homogeneous region are well-known. Therefore, knowing the Interior 
fields at the boundary surface provides the necessary information to solve 
the exterior fields by using the boundary conditions at the boundary surface. 
The finite-elements technique Is used to find the interior field. 

Advantages: 

Adaptable to scatterers of non spherical shapes end Inhomo- 
geneous material. 

Results are baaed on exact Haxwell's equations. 

Disadvantages: 

Requires extensive computer memory. 

Complicated computer codes. 

Based on this assessment It Is noted that only two techniques are 
sufficiently reliable to provide scattering results that possess ilttle 
uncertainty for Pruppacher and Fitter raindrops. They are the extended 
boundary condition method and the finite elesient technique. The complexity 
of the finite element formulation for truly three dimensional problems and 
Its enormous demand for computer memory renders this technique less than 
attractive for extensive computations. The extended boundary condition 
technique Is therefore preferred. 
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III. A Brief Suamary of the Ex::ended Boundary Condition Method . 

The geometry of the problem is shown In Fig. 1. For a linearly polarized 
incident plane wave 

B^(z) - E^e^*“e^ U) 

where is the amplitude of the plane wave, e^ la a unit vector Indicating 
the polarization direction of the Incident wave which la taken to be either 
the X or y direction, k la the free-«pcce wave number, and the wave la 
propagating in the •♦•z direction, the scattered field takes on the following 


form: 

IkP 

E^(r) - f (o, 1) ^ (2) 

where ^ and o ate, respectively, the Incident direction and the observation 
direction. Following the matrix notation of Ishimaru [1978] and Van de Uulst 
[1957], the Incident field and the far-zone scattered field can be written as 
follows : 


li “ =11 ill =1// il// 


( 3 ) 


h " =sx ^1 =s// iis// 

where e^^ or e^^ and e^^ or e^y,^ are, respectively the unit vectors perpendi- 
cular and parallel to the plane of scattering which Is the plane containing 
the vectors ^ and £. Note that is necessarily equal to ^ , but e^yr^ Is 
not necessarily equal to e^// and E^and are related to and E^// by 
the following relation: 
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Tha scattarlng function* 82* and 8^ u«ad by Van da Rulat aa wall aa 
tha 8tokaa paranatara can all ba obtatnad fron f^^* ^]^2* ^21 ^22* 

exanplat 


■ X ^11* “ 7 ^99f ■ 7 ^ 


1 ‘22’ 


c „ JL # 

1 ‘21* “4 1 ^12 


(5) 


and if the incident wave haa an arbitrary atate of polar iaat ion and ita 
8tokea paranatara are given by l2i» 8tokea para- 

neter* l2^, U^, and of tha acatterMl wave are given by 






Re(fl2^fl2) 

-I»<«ll«l2> 

hi 

^2a 

1 


|'22p 

Re(f2if;2> 

-In(f2if22) 

hi 

«s 

■? 

2Re(filf2l) 

2Re(f^2^22^ 

Re(fiif22+f 12^21^ 

*Im(fllf22‘*^12*21^ 

h 



2Im(fjjf2^) 

* 

* * . 

* * 




2Imlf,2»22> 

In(fiif22*t^i2f2i) 

Re(fiif22**^12*2l) 

m 

h 

m m 


( 6 ) 


where * means the compile conjugate. 

Other comtonly uaed parametera nay alao be defined in terms of the 
function f. (o, O; auch aa: 

The differential scattering cross section, 0^(0, O ■ |f <0, 0| 

2 

The bistatic radar cross section, Ohi (o, O ■ 4ir lf(o, pj 

A 

The backscatterlng cross sectlan, 0|, ■ 4ir |f(-^, 0| 


Tlie scattering cross section, 0^ ■ ^^|f (o, _! ) j 
(dO “ differential solid angle) 

The extinction cross section, 0. * o 

A 9 • 


d(i 


(7) 

( 8 ) 
( 9 ) 

( 10 ) 


^ ii (f (it !,))*•£ ■ In(fj^j^(o)) 


■ ^ In(f22(o)) 


( 11 ) 


(0 


^ absorption cross section, e^ ■ unit vector in the direction of 
polarization of the incident wave.) 
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So, knowing ths antrix elemmta fj^2> ^21* ^22’ cnlculat* 

all the dealrad parametera. Theae matrix elementa are obtained by aolvlng the 
boundary valu-i problem. We intend to use the extended boundary condition 
technique to solve this scattering problem. 

The scattering bddy, assumed to be homogeneous and isotropic, is 
characteriztxl by the constitrtive parameters ami t where vi^ is the free space 
permeability and c is the permittivity of the dielectric material (c may be 
complex to account for the lossy case) . The surrounding medium is considered 
to be free space with parameters p^, e^. For the given Incident field, the 
scattered field must be determined. 

The theoretical development for the extended boundary conditon 
technique is as follows: 

First, the Equivalence Principle is applied, breaking the scattering 
problem into two separate problems, an exterior part and an interior part. 

One equation for the unknown field quantities is derived fro^ the external 
problem where it is found that the scattered field due to the surface currents 
must completely cancel the Licident field throughout the Interior volume. 

The internal problem is then considered and the fields within the 
dielectric region are expanded lit regular vector spherical wave functions with 
coefficients to be determined. Superposition is applied and the boundary 
conditions at the surface lead to a linear system of Integral equations for 
the coefficients of the unknown surface fields in terns of the incident field. 

The scattered far field is then determined by evaluating the internal 
field at the surface and substituting into the original expression which gives 
the exterior scattered field in terms of the surface fields. 
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To avoid writing many lengthy mathematical expressions, detailed deri- 
vation will not be given here. We shall simply state the final results: 


(K + */Tj.J) c^ + (L + 1) dy - -Ja^ 

(12) 

V " 1,2, ... ,K 


(I + »^j.L) Cj^ + (J + •'i^K) d^ - -Jb^ 

(13) 


where 

I - — ^ n.H^ (kr') x (k'r')dS 
X S V - ~\i 

J - — ^ n.M^ (kr*) x (k'r')dS 

IT s S> ~ ~\l ~ 

K * — ^ n.N^ (kr’) x (k'r’)dS 
” S V y 

L « — ^ n.N^ (kr’) x (k’r’)dS 
S \J y 

and e “ — . n is a unit vector normal to the boundary surface S and r’ is 
*■0 

the radial vector from the origin to the boundary surface, and are 
expansion coefficients for the Incident wave defined as follows: 

where v is a combined index incorporating o, m and n. is a normalization 
constant 

( "n^l) (nHp) I “3 

” ^m ^(n+1) (n-hn) ! * ® 2 * ^ > 0 

The spherical wave harmonics M and N have been defined by Stratton^^. These 
functions, solutions of the vector wave equation, are given by, 

M^’^(r) - V X r ™^p“(cos 0)z^*^(kr) (15) 

omn sin m* ® 

^’^(r) - i V X »^’^(r) 
omn 1c c^mi 
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where o “ even or odd, 

p"(co8 6) > associated Legendre function, 

1 3 

and 2 ^* (kr) “ an appropriate spherical Bessel function. 

For solutions of the wave equation which must be finite at r ■ 0, s^(hr) 

“ and the resulting vector spherical functions M and N are known 

as solutions of the first kind. The other vector spherical functions we 
will be using are obtained by using the spherical Hankel functions, l.e., 

3 ( 1 ) 

*'^(kr) • h ' (kr) - j (kr) + jn^(kr) . This solution, representing outgoing 
11 “ n u 

waves. Is called the solution of ti.e third kind. 

c^ and d^ are the expansion coefficients for the field Inside the 
dielectric which are defined as follows: 

E(k*r) - Z lcj^(k'r) + d, N^(k’r)] (16) 

jj-1 w-u - w -u - 

where y Incorporates the indices o, m, n and c^ and d^ are unknown 

coefficients, k' ■ and k • ui/y e . The H field Internal to S Is 

r 0 0 

given by 

H(k’r) - -j-i- (7 X E(k’r)) 

- — juy - — 

o 

- ", I%(’ « « Nj<k'r))I 

o y*i 

- - J E (C r(k'r) +d M^(k'r)l (17) 

The coefficients of the scattered fields which are defined as follows: 

E (kr) - E [p M^(kr) + q (kr)] (18) 

V-1 '’V V V 
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tfhere r Is outside a circumscribed sphere, are 


N 

Py E {[K’ + /TyJ’l c + (L’+ ✓Tl’l d ) 

R 

% — JDy Mtl‘ + *^,1'’! Cy + tJ’+ /^K*] d^} 

where 


I' - 4" ^ n.M^(kr’) X M^Ck’r^lS 

J' • 4- ^ n.M^(kr') X N^(kV)dS 
r 8 V ~ -M “ 

K' - — ^ n.N^(kr') x M^(k'r')dS 

IT 8 “ ~\) “ p — 

L’ • V ^ n.N^(kr’) x N^(k'r')dS. 

It 8 — — u — 


( 19 ) 


The exact theoretical problem requlrea the solution of an Infinite 
set of equations. The equations resulting from the External Problem are 
such a set. The first concession to practical considerations Is made In 
the derivation of the Internal Problem, where the Infinite term spherical 
expansion for the Internal field Is truncated to N terms. This finite 
Internal field expansion is then substituted into the first 2N equations 
of the Infinite set of Eqs. (13) and the final set of truncated equations 
is obtained. The scattering results are obtained numerically by solving 
the complete set of equations for successive values of N until the final 
results (such as the differential scattering cross section) converge to a 
specified accuracy. This Insures that enough of the expansion terms have 
been kept to guarantee the correct final result. 


11 



The Indices u and v each run from 1 to N and are related to the 
internal field expansions respectively. As Indicated previously* they 
incorporate the indices o* m and n where a is even or odd and m and n are 
indices of Bessel and Legendre functions. To distinguish between expansions 
over V and y, the indices which they represent are unprimed and primed, 
respectively, i.e., 

V is over o, m, n 

ji is over o' , m* , n' 

where the expansion scheme is as follows: 

omn " eOl, oOI, ell, oil, e02, o02, el2, ol2, e22, o22, ... 

where e “ even and o “ odd. 

There are essentially four main steps in the solution for the diffe- 
rential scattering cross section. 

(1) Solution of the system of Eqs. (13) for the coefficients of the 
internal field c^ and d^ in terms of the Incident field coefficients a^ 
and and S is the surface of the scattering object. Note that the 
solution of the Eqs. (13) requires a matrix inversion procedure. 

(2) Substitution of the coefficients c^ and d^ into Eqs. (19) to 
determine the coefficients of the scattered field. 

(3) Substitution of the scattered field coefficients into a spherical 
expansion for the far field vector amplitude. 

“it ^ [p ^ (e.4»)+ Jq B (e.*) 

“‘’11 ja \T^ V V 

(A) Computation of the differential scattering cross section 
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IV. Display of Results and Discussion 

Selected Pruppacher-Pltter raindrop shapes are displayed in Fig. 2. 

It can be seen that as the raindrop size Increases, Its shape deviates 
elgnlflcantly from a sphere and It even differs noticeably from an oblate 
spheroid . These Pruppacher-Pltter raindrop shapes were used in our compu- 
tation and compilation of the scattering results. Specifically, we have 
obtained the scattering matrix elements fj^ 2 > ^21* ^22 following 

parameters : 

Raindrop sizes: 0.25 mm (0.25 mm) 3.5 mm 

Frequency: 30 CSz 

Temperature: 20*^C 

e : o'- (10°) 360° 

8 

a : 0° (30*> 90° 

^ ; 00 (30°) 180° 
s 

Results are tabulated in Table 1 which has been condensed into the two 

microfiches attached at the end of this report. 

We have, of course, compared our results with those obtained by Cross 
4 5 

and Morrison , and Oguchl who used the point-matching technique and by 
12 

Morgan who used the finite-element approach. The results all agree within 
a few percentage points. These differences may be caused by numerical 
errors due to finite grid size, finite number of matching points, finite 
matrix size, finite integration steps, accumulated numerical errors or the 
limitation of the particular numerical scheme. Based on this comparison, 
one may conclude that the point -matching technique actually provides rather 
good results in spite of uncertain convergence and the lack of rigor due to 
"Rayleigh Phenomenon." It is also significant to note that our results agree 
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quite well with those obtained according to the finite-element approach which 
Is based on the exact Maxwell equations. 

It Is of Interest to note that the scattering results for a Pruppacher- 
Pltter shaped raindrop compared very closely to those for an equivalent 
oblate spheroidal shaped raindrop even when their shapes are visibly different 
for large size drops (see Fig. 2). It appears that the absorptive character- 
istics of water at 30 GHz and 20^ C are responsible for minimizing the 
differences In the scattering results due to shape differences in the Pruppacher- 
Pltter and equivalent oblate spheroidal drops. 

To illustrate the effect of drop size changes on the scattered field. 

Figs. 3-5 are introduced. In Fig. 3 the normalized scattering, absorption and 
extinction cross-sections are plotted as a function of raindrop sizes. The 
relatively smooth character of the curves can again be attributed to the 
absorptive nature of water. Polar plots of scattered fields In the * 90° 

plane for canting angles o ■ 0° and 63° and for various drop sizes are shown 
in Figs. 4-6. It can be seen that no cross polarized field is generated 
for 8ir<all drop sizes (with spherical shapes) or for canting angle a ■ 0. 

(See Fig. 4.) For large drop sizes (with non spherical shapes) and for 
canting angle o ^ 0, a cross polarized field is generated. (See Fig. b and 6.) 
One also notes that the induced cross polarized field Is at least one order 
of magnitude weaker than the non-cross polarized scattered field. This is 
again due to the absorptive characteristics of the raindrops under consideration. 

In conclusion one may state that the extended boundary condition method 
Is an efficient method in yielding reliable scattering results for non- 
spherlcal Pruppacher -Fitter raindrops. One notes that the large absorption 
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constant of water at 30 GHz and 20° C tends to minimize the shape effect 
on the scattered field. It is expected that, due to the transparent nature 
of ice at 30 GHz, the shape factor for ice particles will play a significant 
role in giving rise to the generation of cross polarized field. 
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a « CANTING ANGLE 




Figure 1. The geometry of the problem. 
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PARAMETERS FOR CALCULATIONS 

PRUPPACHER • PIHER RAINDROP SIZES 0. 25 mm (0. 25 mm) 3. 5 mm 

FREQUENCY: 30 GHz 
TEMPERATURE: 2fc 
a = (f (3(f)90l® 

= 0® (30®) 180® 

= 0®(i0®)360® 



0.25 mm 




RAINDROP SHAPES 

Figure 2. Pruppacher-Pltter raindrops 
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^40RMALIZED CROSS SECTION 



Figure 3. Normalized cross section as a function of raindrop size 
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Ift CM* 


Figure 6r Polar dlagraas for scattered fields with a ” 60” (cross 

polarization present for non-spherical raindrops). Incident 
field is polarized in the x-direction. 



